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Abstract 

Testing of the gravitation equations, proposed by one of tlie authors earlier, by a 
binary pulsar is considered. It has been shown that the formulas for the gravitation 
radiation of the system resulting from the equations do not contradict the available 
observations data 

1 Introduction 

In paper [|l| gravitation equations which do not lead to a physical singularity in the 
center of the spherically symmetric field were proposed. These equations also predict 
that there can exist stable supermassive compact configurations of the degeneration 
Fermi - gas without an events horizon The equations do not contradict classical 
tests at the distances from the center which are much larger than the Schwarzschild 
radius. In the present paper we find the power of the gravitational-wave radiation from 
a close binary system and use the result to find the deceleration of the orbital period of 
the pulsar PSR1913+16 conditioned by the gravitational radiation. In this case we deal 
with a moderately strong gravitation field and use the definition of gravitational energy 
that follows from the gravitation equations under consideration. 

2 Equations of Field 

Thirring [||] proposed that gravitation can be described as a tensor field ipafi {x) of spin 
two in Pseudo-Euclidean space-time E4 where the Lagrangian action describing the 
motion of test particles in a given field is of the form 

L = -mpc[gapirl^)x°'x'^Y^'^ . (1) 

In this equation gap is a tensor function of i/jap, rrip is the mass of the particle, c is the 
speed of light and x°' = dx°' /dt. 

A theory based on that action must be invariant under the gauge transformations 
tpa0 — * that are a consequence of the existence of "extra" components of the 
tensor Transformations ^pa/s — * "00/3 give rise to some transformations gap — > 
(jap- Therefore, the field equations for ga/six) and equations of the motion of the test 
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particle must be invariant under these transformations of the tensor gap. Equations 
of gravitation that are invariant with respect to arbitrary gauge transformations were 
proposed in [111] . These equations are of the form 



r>7 

a/3;7 



0. 



The equations are vacuum bimetric equations for the tensor 



(2) 



(Greek indices run from to 3), where 



n. 



Q/3- 



(3) 



Kp-in + l)-' 
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are the Christoffel symbols of the pseudo-Euclidean space-time E4 whose funda- 
mental tensor is rjap, T'^p are the Christoffel symbols of the Riemannian space-time 
V4 with dimension n — 4, whose fundamental tensor is gap. The semi-colon in eqs. 
denotes the co variant differentiation in E^. 

A peculiarity of eqs.(^ is that they are invariant under arbitrary transformations 
of the tensor gap retaining invariant the equations of motion of a test particle, i.e. 
geodesies in V4. In other words, the equations are geodesic-invariant. Thus, the tensor 
field gap is defined up to geodesic mappings of space-time V4 (in the analogous way 
as the potential Aa in electrodynamics is determined up to gauge transformations). A 
physical sense have only geodesic invariant values. The simplest object of that kind is 
the object B^p which can be named the strength tensor of the gravitation field. The 
coordinate system is defined by the used measurement instruments and is given. 

The Christoffel symbols are transformed under the geodetic (i.e. projective) map- 
pings as follows: 

f;;^ - + ^p6^ + ip^s'^p, (6) 

where tp^ is a vector-function of x". This equation has a simple interpretation in an 
(nH-l)-dimensional manifold M.^ where the admissible coordinates transformations are 
of the form 

-a -a, 1 2 3\ rn\ 

x=x{x,x,x,x), (7) 



-4 4 

X = X 



(fadx" 



(8) 



Namely, eq. (^ can be interpreted as the transformation of 4-components T^^ of the 
connection coefficient (A, B,C 0..4) in M5 under the transformation (^ if 
the condition T^p = is satisfied. 

For this reason we will consider the tensor gap as 4-components of 5-dimensional 
tensor 

gap gai 

54a 544 



gAB 



(9) 



The components are transformed under (^ as follows 

gap = gap + ga'-Pa'^P + gia'-PP + g4p'Pa, 



(10) 
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9aA= 9ai+ giWa, (H) 
.944 =344- (12) 

Eqs. ( [lT| ) coincides with the transformation of the object (or Qq, = r(^^ — r^,^) 
under (g), if 544 = n + 1. For this reason, we will assume that 

We will assume also that ( p^ is the transformation of the tensor gap under the 
geodesic mappings of V4. Then there exists the geodesic-invariant tensor 

Ga/B = - (u + ly^QaQp, (14) 

and the geodesic-invariant generalization of the Einstein equations with matter source 
are of the form 

- BZaBfj^ = k {T^p - l/IG^pT) , (15) 

where k — SttG/c^, T^p is the matter energy-momentum tensor, T — G'^'^Tap. At 
the gauge conditions Qa = they coincide with the Einstein equations. 

Consider now the question about the definition of the energy of gravitational field 
in the used theory. Set in eqs. ( |l5| ) Tap = pc^UaUp, where p is the matter density and 
Ua is the 4-velocity of matter points. At the small macroscopic velocities of the matter 
we can set uq = 1 and Ui = 0. Then, the 00-component of eq ( [1^ can be written in 
the form 

^oo;/3 = x(pc'+ioo), (16) 
where x — k/2, and too is the 00-component of the tensor 



Setting 

we obtain the equation 



= x-'Bl.B'^p^. (17) 
= c-^dU/dx" (18) 
AU = xipc^ +too)- (19) 



In the absence of the term too and at the distances from the central masses much larger 
than the Schwarzschild radius this equation coincides with the Poisson equation for the 
Newtonian gravity potential. For this reason it is natural to expect that in general case 
eq. ( [l9| ) is the differential equation for a generalization of the gravitational potential 
too and the term too, which is in the equation additively with pc^, is the energy of the 
gravitational field. 

To verify this assumption let us find the energy of gravitational field of the point 
mass M as the following integral in the Pseudo-Euclidean space-time: 

S = J toodV. (20) 

In the Newtonian theory this integral is divergent. In our case we have 

^00 = 2x^^Sdo ^01 (21) 
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and, therefore, using spherical coordinates, we find 



where 



and 



1 r'^c^ 

toodV = -^J, (22) 



B(z, w) = / ^—dt (24) 



is B-function. Using the equaUty 

?l 7- in) = 

T{z + w) 



(25) 



where T is F-function we obtain 
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£ = ^ = Mc^ (26) 

We arrive at the conclusion that the energy of the point mass is finite and the rest 
energy of the point particle is caused by its gravitational field. The spacial components 
of the vector Pa — toa are equal to zero. Due to these facts we may consider (at 
least in a weak-field approximation) the tensor <a/3 as the energy-momentum tensor of 
gravitational field. 



3 Gravitation radiation of a binary system 

In no gravitation case g^p = rj^fj, where ?7c,/3 is the metric tensor of Pseudo-Euclidean 
space-time £^4. For this reason it is natural to suppose that in a weak- field approxima- 
tion 

gaf} = Va/S + hap. (27) 

To simplify eq. (^Sj) let us choose the following gauge condition 

r^'^PDip+'^r^D^^^^Q, (28) 
^ n + 1 

where D'^^ is the Christoffel symbols of space-time V4, in which the derivatives are 
replaced by the covariant ones in i?4. Equality (GS) is covariant and, therefore, does 



not imply restrictions for choosing the coordinate system. Using (27) we obtain to the 
first order in hajs: 

Blp = I {hip + hlo. - h'^p) - ^ (" + 1)"' {S1K;P + S;hQ ■ (29) 
and the gauge conditions (^8|): 

hip = {n+iy'K.,^. (30) 

Consider a system of the slowly moving bodies of a finite volume where T^^, 7^ 
in the vicinity of the coordinates origin, and suppose that jTool ^ \TQi\ S> \Tij \ (Latin 
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indexes run from 1 to 3). Using an orthogonal system of coordinates with the metric 
Vfiv = diag{—l, 1, 1, 1) let us find the gravitation radiation in the wave zone. 
Equations ( p^ ) now are of the form 

□V^«/3 = -2kTo,0, (31) 

where 

= hai3 ~ l/2riaph\. (32) 
The general solution of eq. ( ^ ) for departing waves is of the form 

. f r„,(c, I. ..I.r-l ^^.^,^ ,33) 

2tt J |r-r'| 

where r ' is the radius-vector of an infinitely small element of the source and r is the 
radius vector of the point in which the field is calculated. 

If the origin of the coordinate system coincides with the masses center of the radi- 
ating system, the solution of eq. ( ^1| ) is given by 

= const + o(r"^), = o(r"^) (34) 

^P^j = -^i{ct - r, x^) + o(r-2), (35) 



where 



hj ^ ^ J Toox^x'^d^^'^x', (36) 

the points denote the differentiation with respect to time t and const is denotes an 
additional term to the newtonian potential. Therefore, the solution of eq. ( ^3| ) for weak 
gravitational waves can be written in the form 

V'o,. - 0, (37) 



^zfc = T^4(ci-r,x"=). (38) 

Therefore, the lowest mode of the gravitational radiation is the quadrupole one. 
Let us find the number of the independent states of the field ijjap in vacuum. Ac- 
cording to ( ^I| ) the functions ^af3 in vacuum obey the equation 

□^c./3 = (39) 

and gauge conditions (^0|), which give four constraint equations. Equations ( ^7| ) also 
gives four constraints. Taking also into account that for weak waves tpoi = 0, we have 
only two independent components of ipap- It follows also from ( ^0| ) and (^7|) that the 
trace of the tensor ipap is equal to zero. Thus, 

^om = 0> -^tj , J = 0, V'fcfc = 0, (40) 

where the coma denotes the partial derivative. Therefore, the plane gravitational wave 
is a transversal with the zero trace. 

Now let us find the flow of energy of the gravitational radiation. Starting from 
eq.(p^, we find that the flow of energy carried by the wave in the radial direction is 
given by 

-ctok = -rri'ipiji^tj), (41) 
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where n'^ — /r and the symbols ( ) denote the period T-averaging of the a function 
f{t) by the formula 

(/W> = ^^ f{t)dt. (42) 

We can also satisfy condition ( |40| ) by projecting the tensor t})ap to the plane perpendic- 
ular to the vector with the help of the operator 

Pij — 5ij — niTij. (43) 

As a result, the transversal-traceless part of the tensor ipap is 

i^I^^ = ^rPjsAs - \P^J {PrsAs) , (44) 

It should also be noted that the transversal-traceless part of the quadrupole momen- 
tum lij, in (^^, and the magnitude = — l/3SijIkk are the same, which can be 
verified immediately. For this reason 

= (45) 

Finally, we obtain 

rirk 1 

~'^^0r = g^^2^2g ^'^'3 '^»j ~ 2JikJkjnin.j + -JtjJkiniTijnuni), (46) 

The averaging time rate of the energy flow change can be found by integration of 
eq.(^) over the sphere of the radius r: 

= / ctor-rirr'^ sin 9 dO dtp = -^(Jijhj), (47) 

were E is the total mechanical energy of the radiating system. This formula completely 
coincides with the one in general relativity. 

Consider a system of two gravitationally bound point masses mi and m2- As- 
suming that the distance between the bodies is much larger than the Schwarzschild 
radiuses we can use the Newtonian mechanics. Since E = — G/iAf/(2a), where 
M — mi + 7712, M = mim2/M, a is the large semiaxes of the orbit and the or- 
bital period is Pi, = 27ra'^/^(M G)~^/^, the relative velocity of increasing of the orbital 
period is given by 

dPb/dt _3da/dt _ 3 dE/dt 



Pb 2 a 2 E 

It gives 

dPb/dt 1927r 



(48) 



Pb 5c5 



(27rG)^/^ F^^/^ /(e)mim2Af-^/3 , (49) 



where /(e) = (l - e^) (l + 73/24 + 37/96 e''). It should be noted that in the 
above calculations the shift of the periastron is not taken into account since at the used 
one period-averaged the shift is very small. 
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4 Application to PSR1913+16 



Let us use the above results for an analysis of the binary pulsar PSR 1913+16 in the 
framework of the Blandford-Teukolsky model |Q|. Namely, we take into account the 
dependence of the phase of the detected radiation on the following observable param- 
eters that are the functions of the unknown masses mi and TO2. 

1. Increasing the orbital period because of the gravitational radiation ( p9|) . 

2. The rate of the periastron shift of the pulsar orbit which according to [H| is given 

by 

. Q^GM nG^AP 

w = — ; H T 11 (e), (50) 

where /i (e) = (54 + 16e^ — e^) . The second term for the considered system is about 
0.01% of the total value of uj and, therefore must be taken into account. 

3. The magnitude 

7=—^, (51) 
sin(£;) 

where ti is the moment of the radiation of the pulse, measured in the inertial frame of 
reference of the distant observer, t2 is the same moment, measured in the proper frame 
of reference of the pulsar, E is the eccentric anomaly. 

Supposing that 4-vector k" — {cj/c, k} of a photon, the motion of which is de- 
scribed by the Lagrangian (|l]), satisfies the equality gc/sfc^fc'' = 0, we obtain approxi- 
mately fc" = cj/ Cy/gao, where uj is the frequency of the electromagnetic wave and 500 
is given by [|l]] 

goo = l~rg/irl+r^y/^ (52) 

Let Vp be the pulsar velocity relative to the masses center of the system. Then we 
obtain the following relation between the measured frequencies of the signal cuq and uj 
in the frames of the observer and pulsar 

1 - jvp/c) cos(t^) 
Wo , = — , (JJ) 



VI - rg/fc 1 1 - vl/c 



where d is the angle between the wave direction and the direction of the motion of 
the source radiation, fc = {r^ + VgY^'^ , Tg is the gravitational radius of the pulsar 
companion, r is the distance between the pulsar and its companion. From ( p3| ) we 
obtain, setting -d ~ 7r/2, the relation between time intervals dr and dt in the reference 
frames of the pulsar and observer 

which with the used accuracy coincides with the same relation in general relativity. 

The Shapiro effect is too small (the same as in general relativity) and is not taken 
into account in the model under consideration. 

Proceeding from these results, we obtain in an analogy with [Q] 

GPhern^im^ + 2too) 
7 = — — — —. (55) 

27rac2(mi + 7712)2 

Fig. 1 shows that in the plane TO1TO2 the curves Pb = Pb (rni. 771,2), w = w(m,i.TO2), 
7 = 7(7711^7(1,2) intersect in one point that means that the theory does not contradict the 
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observation data. The value of the full system mass M, resulting from the measured 
value of Pb and e by eq.@ is equal to (2.82845 ± O.OOOO4)M0 and the masses of 
the pulsar and its companion, found by (^5]), are equal to (1.441 ± O.OO3)Af0 and 
(1.387 ± O.OO3)Af0, respectively. These results differ very little from the ones result- 
ing from general relativity gl ((1.442 ±O.OO3)A/0 and (1.386 ±0.003)Mq). Due to a 
kinematic effect in our Galaxy [|] the small correction (-0.017 ± 0.005) • lO'^^ 

must 

be added to the found theoretical value of Pfo = (-2.40249 ± 0.00029) • lO^^^ Taking 

into account this correction we obtain that the ratio of the observational value of Pb to 
the found theoretically is equal to 1.0023 ± 0.0047. 




1 m1 2 3 



Fig.l The plots of the functions 7 = 7(7711,7712) (curve 1), u — w(77ii, 7112) (curve 
2), Pb = P6(mi, 7712) (curve 3). 
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